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Abstract 

We study the interactions of a new spin-1 mediator that connects the Stan- 
dard Model to dark matter. We constrain its decay channels using monojet and 
monophoton searches, as well as searches for resonances in dijet, dilepton and 
diboson final states including those involving a possible Higgs. We then interpret 
the resulting limits as bounds on the cross-section for dark matter direct detec- 
tion without the need to specify a particular model. For mediator masses between 
300 and 1000 GeV these bounds are considerably stronger than the ones obtained 
under the assumption that the mediator can be integrated out. 
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1 Introduction 



Although the Large Hadron Collider (LHC) has yet to see any evidence for physics 
beyond the Standard Model (SM), we know from astrophysical and cosmological obser- 
vations that the SM is incomplete, because it lacks an adequate candidate for the dark 
matter (DM) particle. Moreover, there are likely to be new interactions that connect 
the DM particle to the SM. In this paper, we discuss how colliders - and especially the 
LHC - can contribute to the search for DM and such new interactions. 

The search for dark matter has mainly been focused on particles with mass and 
interactions set by the Fermi scale, e.g. neutralinos in extensions of the SM with softly 
broken supersymmetry (SUSY) at this scale. With an impressive increase of sensitivity 
over the past few years, direct detection experiments such as XENON100 |1| and CDMS- 
II [2j|3) have pushed down upper bounds on the (spin-independent) scattering cross- 
section on nuclei for such particles and have begun to constrain the relevant SUSY 
parameter space. 

For much lighter DM particles, however, direct detection bounds become significantly 
weaker and at the same time such particles must have stronger interactions if their relic 
annihilations are to result in an acceptable DM abundance. Hence an improvement of the 
bounds for light DM is of great interest especially since recent results from the DAMA |1] , 
CoGeNT [5] and CRESST-II |6] experiments hint at DM with mass m x ~ 10 GeV and 
cross-section a p ~ 10~ 40 cm 2 . If DM does have such properties, then production of DM 
pairs at the LHC would be sizeable and result in a variety of observable signals. 

Among the most promising signatures of DM at high energy colliders are excesses 
of events with either a single high-energy jet or a single high-energy photon and a large 
amount of missing transverse energy (MET). Such monojet and monophoton searches 
have been performed at LEP, Tevatron and the LHC but no excess has been observed 
over expected SM backgrounds [7 12 . If the mediator of the DM interaction with the SM 
is so heavy that it cannot be produced on-shell at the LHC, then these searches directly 
bound the coupling of DM to nucleons and are competitive with the bounds set by direct 



detection experiments 13-21 . However, for such large mediator masses, the relevant 



cross-sections for both the LHC and direct detection experiments will be very small 



unless the coupling constants approach the bounds from perturbative unitarity 18 -20 



Of course the mediator mass may well be comparable to LHC energies. In this case 
an effective operator description is no longer valid, because the LHC can resolve the 
interaction and produce the mediator on-shell, which complicates the comparison with 
direct detection experiments. On the plus side, it opens up the possibility to search for 
resonances in various channels from the decays of the new mediator into SM particles. 
Combining the limits from all relevant collider searches, it is still possible to constrain 
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the direct detection cross-section in a model independent way. 

In particular, if DM is light compared to the dominant mediator, we can bound 
the DM direct detection cross-section in terms of the total width and the invisible 
branching ratio of the mediator in a simple way. These two quantities may in turn be 
constrained by collider searches without having to specify an underlying model - even 
though the collider bounds may be much stronger in a specific model framework. In 
this paper we apply this approach, assuming that the interaction between DM and the 
SM is dominated by the exchange of a neutral spin-1 state, here termed R. An example 



of such a spin-1 state is the Z' associated with a new broken U(l) symmetry 22 -35 



Another example is a new resonance associated with a strongly interacting extension 



of the SM 36 37 , e.g. an analogue of the neutral isospin zero u> resonance in QCD. 



Recently, there have been analyses of the 'dark Higgs' associated with the Z' 38 and 



LHC signatures of a baryonic Z' 39 



The outline of this paper is as follows: In Section [2] we introduce an effective La- 
grangian for the spin-1 state, discuss the decay channels and present the production 
cross-sections at LHC and the Tevatron. In Section [3] we first summarize the resulting 
collider bounds from LHC and Tevatron on the various decay modes of R, before dis- 
cussing the bound for each decay mode in more detail. In Section [4j we compare our 
results to limits from direct detection on spin-independent and spin-dependent inter- 
actions. We do this both model independently and in the framework of a spin-1 state 
coupling to the SM via kinetic and mass mixing only, as discussed above. A discussion 
of all assumptions and their validity is given in Section [5] together with our conclusions. 
Appendices [A] and [B] provide, respectively, all relevant formulae for the partial decay 
widths and the coupling constants of R. 



2 Interactions of a neutral spin-1 mediator 

We start from an effective Lagrangian (similar to [28] ) describing the interactions of 
the neutral spin-1 state R with the SM fields and the DM particle. We then discuss all 
possible decay channels and describe how they can be constrained by collider searches. 

2.1 Effective Lagrangian description 

We divide the Lagrangian into the couplings to DM, SM fermions, SM gauge bosons, 
the Higgs, and anything else 

C R = + £// + £gaugc + £-h + £-x ■ (1) 
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Depending on whether DM is a Dirac fermion \ or complex (pseudo-) scalar 0, we define 
= or = as appropriate, where 

£* = R»xY(9^R - ^7 5 )X , Cf = g^Jf (2) 

and J% = i((j)*d^(f)- (pd^cj)*). We will not consider the CP-odd operator P^<9 M (0*0) here 
(for a discussion see e.g. [40]). 

The interactions of i? with the SM fermions are described by 

C ff = E ^/V(& - ^ 7 5 )/ , (3) 

where q, £, v denote SM quarks, charged leptons and neutrinos respectively. 

Neglecting CP-violating terms (see e.g. 41 for a more complete discussion) the 
couplings of R to SM gauge fields can be written as 

£?auge = g^wAiRW+w-^ + g^Mw+w-jh 
+ g% wwl ((RZW + w-)) + £ wwi ((R'yW + w-)) 

+ g B zz [[RZZ}} t + g%\\RZi\\ t + g% W3 [[RW + W-]] e 

+ g§ WW2 e^R,Z u w;w; + g* WW2 e» vp ° 'R^W+W; , (4) 



where 



[[RW+W-]}! = i[(d IM W+-d v W+)W f "-R v -(d ll W--d u W-)W fi+ R v ] 
[[RW + W~}] 2 = faRv - dvRjiW+W- -W»~W U+ ) , 



[[RV x V 2 ]] e = e^(V lfl d p V 2v -d p V lfM V 2v )R 
((RVW + W-)) = 2R i y^W~W u+ - R^W^ + V V W U - - R^W^V U W U+ , 

for appropriate combinations of Vi — {7, Z,W + ,W~}. The operators in the first two 
lines of Equation (HI) conserve C and P separately, while the operators in the last two 
lines are CP even but P odd. If the underlying theory violates CP, then CP violating 



couplings are also possible 42 . In principle there could also be a coupling between 
the R and two photons, leading to the decay R — > 7*7. 1 Such a decay would however 
correspond to a non-renormalisable operator which - in order to have any observable 
effects - would imply additional new physics at rather low scales. We will not discuss 
such operators. 



1 Note that the Landau- Yang theorem only applies to on-shell photons. 
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For comparable couplings the triboson final states are suppressed compared to the 
diboson ones due to smaller available phase-space This is expected when R is a new 
gauge boson (for details see e.g. (43)). Consequently, we will neglect triboson decays of 
R in the following. 2 

Finally, couplings of R to the SM Higgs are of the form 

£h = g% H R»Z»H + g§ HH R,Z»H 2 . (5) 

Again, we expect decays into ZHH to be significantly suppressed compared to decays 
into ZH so we neglect them in the following. A coupling of R to HH is absent because 
if is a real scalar and terms proportional to d^W 1 are CP violating hence we neglect 



them 41 



Leaving L\ unspecified for now, the decay modes of the vector R may then be 
summarized as 

r fl = r** + ^r 99 " + ^r" + ^r^ + r w + r zz + r 7Z + r ZH + r x , (6) 

q l v 

where the formulae for the partial widths are provided in Appendix [Aj For consistency of 
our description, we will impose V r/itlr < 1 which already gives a bound on all coupling 

constants. If we define = J {g^ R ) 2 + (g$n) 2 for any fermion ip we then have the 
following constraints on the couplings in isolation 

g<f>R ^ 12 , g x < 6 , g e , g v < 3.5 , g q < 1.5 . (7) 

Here we assumed family- in dependent SM couplings: g u = gj = g q . 



2.2 Production and decay of R at colliders 

At colliders, the new spin-1 state R can e.g. be produced via Drell-Yan (DY) production, 
vector boson fusion (VBF) or 'i?-Strahlung' from a SM gauge boson. We will focus in 
this paper on DY production. However, VBF will be important if the coupling of R to 
W's is large, e.g. if R arises from a composite theory. In this case, there would be more 
search channels with 2 additional jets in the final state. 

In the case of DY production, we can decompose the cross-section for the production 
of R in association with an additional particle Y and subsequent decay of R into xy, as 

a(qq ^R + Y^xy + Y)= a(qq -> R + Y)- BR(R -> xy) , (8) 

2 However if R is a 'techni-omega', the coupling to the triboson final states can be enhanced so that 
the corresponding decays contribute significantly to the total width of R. 
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Figure 1: Drell-Yan production of R as a function of rriR at the LHC (with y/s = 7 TeV) 

and at the Tevatron for equal to those of the SM Z boson (solid lines). We also show 
the production cross-section for the case that R couples either only to u-quarks or only to 
d-quarks. 



where we have used the narrow width approximation (NWA), applicable if Th/w,r <C 1. 
As a consequence of the NWA, the DY cross-section a(qq — > R — > xy) can be written 

as 



44,45 



a xy oc [glw u (s, m 2 R ) + g 2 d w d {s, m 2 R )] ■ BR(R -» xy) , (9) 

where w u ^ d parameterise the parton distribution functions (PDFs) of the proton. For 
a narrow resonance the only dependence of the w u , d coefficient on the resonance R is 
through rriR. Since the functions w u>d are known, one can always translate a bound on 
g u into a bound on g d or a bound on a model with a given ratio of g u /g<i into a bound 
on a different model. For this reason, we will show experimental bounds assuming 
9u = 9d = 9q fo r simplicity. 

In Figure [I] we show the leading order (LO) cross-section for a(pp — > R), with 
the couplings of R to SM fields set equal to those of the SM Z boson as well as the 
individual contributions from w-quarks and (i-quarks. The production cross-section of R 
with different couplings may be found by a simple rescaling. Note that QCD corrections 
can enhance the DY production significantly; in the following, we take the K-factors 



from [45] for the MSTW08 NNLO PDF. 

Once we know the production cross-section of R for given couplings g q , we can 
translate LHC bounds on the cross-section for a certain final state xy into a bound on the 
product g 2 -BR(R — » xy). The decay modes of R into the X sector can appear at colliders 
as additional missing energy, displaced vertices or high-multiplicity SM final states. An 
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Figure 2: Bounds on ■ BR(R — > xy) as 
a function of mR, with xy being either SM 
gauge bosons or ZH. 



Figure 3: Bounds on ■ BR(R — > xy) as 
a function of m^, with xy being either SM 
fermions or DM particles. 



example is the decay mode R — > Zh! where h! is a new scalar state responsible for the 
mass of R. Another contribution to T x could come from decay modes of R to additional 
hidden sector states. A number of such possibilities were considered e.g. in 46 



3 Collider bounds 

In this section we summarize the collider bounds from LHC and Tevatron for the de- 
cay modes of R and then discuss each decay mode in more detail. Table [T] lists the 
current collider searches we consider here, while the corresponding limits we derive are 
summarised in Figure [2] and Figure [3] The confidence level for all bounds is at least 
95%. 



3.1 Mono jet searches 

The monojet final state arises from production of R with an additional jet, j = q,g, 
followed by the decay of R into DM or neutrinos: a{pp — > jR — > jft T )- Decay modes 
into additional hidden sector states also contribute to the monojet signal, provided these 
states do not decay back to SM states within the detector. 

Monojet searches have been performed at the Tevatron [8,58 , and at the LHC by 
both CMS [9l and ATLAS |59l 



with similar sensitivity. To calculate our bounds we 
compare the limits from ATLAS with the parton-level monojet signal from R simulated 
using CalcHEP 60 . For monojet searches, the jets have sufficiently high that the 



errors from neglecting parton showering and hadronization are small (see e.g. 14,61]). 
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Channel [Exp] 


t r pi — "1 -i 

L [fb ] 


Mass range 


Couplings 


Reference 


v A i n*k f ATT A Cl 

PP -> 3 + f T [Al LAS] 


l.U 


* 










PP^3+f T [CMS] 


4.7 


* 


9g9 x 




n 

9 




PP ^ 3 +f T [^Ut j 


0. ( 


* 


9q 9x 




8 




v — , i l~/^T\ TCl 

-» 7 + f T [UMSJ 


11/1 

1.14 


* 


9q9x> 9q9z 1 




11 




PP -> 7 + |* T [CMS] 


4.7 


* 


9q9 X i 9q9% 




12 




k A A \ ATT A Cl 

PP ^ 3 3 [A1LASJ 


1 n /i 
1.U4 


nnn /innn 

yuu — 4uuu 


9q 




47 




PP -> J J [CDF] 


1 1 Q 

l.lo 


o c n i/i nn 
ZOO — 14UU 


9q 




48 




-»■ [CMS] 


4.9 


300 - 2500 


9q9e 




49 




~™ v [ATT A Cl 

pp —> it [A I LAb\ 


i no i o 1 
1.08-1.21 


200 — 2000 


9q9t 




50 




pp -> rr [CMS] 


4.9 


350 - 1600 


9q 9r 




51 




pp -> Z Z ( *> [ATLAS] 


4.9 


110 - 600 


9q 9zz 




52 




mm V V 7 [ATT A Cl 

pp ^ Z Z [AlLAoJ 


1 n 
l.U 


oZU — IbUU 


9q 9zz 




53 




pp — > Z Z K ' [OIVISJ 


/l fi /I fi 
4.0 — 4.5 


1 1 n finn 
11U — ouu 


9q 9zz 




54 




pp ->■ W [ATLAS] 


4.9 


110-600 


9q 9ww 




52 




pp^W WW [CMS] 


4.6-4.8 


110-600 


9q 9ww 




54 




pp -»■ Z/T [ATLAS] 


4.7 


* 


9q 9zh 




55 




pp^rtt (boosted) [CMS] 


4.6 


1000 - 3000 


9q 9t 




56 




pp — >• ti (semileptonic) [CMS] 


4.9 


500 - 1500 


9q 9t 




57 





Table 1: Collider searches in final states that constrain the couplings of R. Fields marked with 
a * correspond to searches that do not look for heavy resonances and which, consequently, give 
constraints for arbitrary wir. In the case of scalar DM g x should be replaced by g&R. 
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Figure 4: Limits on ■ BR(i? — > inv) 
from monojet searches (red, dotted and pur- 
ple, dot-dashed) [lOj and the monophoton 
searches (green, dashed) [TJ] . The red dotted 
line corresponds to the cut pt > 350 GeV, 
the purple dot-dashed line to pp > 250 GeV. 



Figure 5: Limits on g^-BR(R — > Z7) from the 
monophoton searches in Ref. [IT] (red, dot- 
ted) and Ref. [TJ] (light blue, solid). 



For m R > 400 GeV, we find that the ATLAS search with prC?) > 350 GeV gives the 
strongest constraint. For lighter R, a stronger bound is obtained from the ATLAS search 
with pt(j) > 250 GeV. The Tevatron gives the strongest bound for uir < 100 GeV be- 
cause of its high luminosity and lower monojet pt cut. Our results are shown in Figure |4| 



3.2 Monophoton searches 

The monophoton final state can arise from two different processes. The first is similar to 
the monojet process with production of R and initial state radiation of a photon, followed 
by decay of R into DM or neutrinos a{pp — > jR —> jf T )- The second possibility is DY 
production of R followed by the direct decay of R into Zj, with subsequent decay of Z 
to neutrinos a{pp — > R — > jZ — > 7 vv). 

We find that the initial state radiation of a photon provides weaker constraints on 
the invisible branching ratio of R than the initial state radiation of a jet (see Figure |1J). 
The second process, however, offers an interesting possibility to limit the branching ratio 
of R — > Zj. Note, however, that such monophoton searches are only sensitive to these 
decays if m#/2 is sufficiently larger than any cut on ft T or the photon p^. We use the 



data from 11,12 and obtain the limit curves shown in Figure |5j using CalcHEP. Note 
that the limits on the direct decay R — > Z'j have been calculated without including the 
contribution of invisible decays of R together with a photon (from initial state radiation). 
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Figure 6: Combined dijet limits from CDF 
(blue, solid) and ATLAS (green, dashed). The 
line width reflects the dependence of the AT- 
LAS bound on T/j (which is varied between 
10% and 25% of the mediator mass). 



Figure 7: Limit on g 2 q ■ BR(i? -> tt) 
from searches for boosted tops (blue, dot- 
dashed) 1 56 1 and semileptonic tops (orange, 
dotted) [57]. 



For sizeable invisible branching of R, the bounds would become even stronger. 
3.3 Dijet resonances 

Searches for resonances in the invariant mass distribution of dijet events have been 



carried out at Tevatron and at the LHC We use the CDF limit 48 for tur < 900 GeV 



and the ATLAS limits 47 for ran > 900 GeV. In the latter case, the bound on the 
cross-section is not quite independent of the width of the resonance, which depends on 
T R and on the detector resolution (see also |61|). 



To estimate this dependence we have generated m,-,- distributions for different val- 



ues of Tr using LanHEP 62 and CalcHEP and convoluted these distributions with the 
detector resolution (47|. By comparing the resulting width of the peak to the bounds 
given in Table II of |47| we estimate how the limit varies with the mediator width for 
Tn/mji in the range 0.1- 0.25. The result is shown in Figure [6] together with the limits 
from dijet searches from CDF. We take the upper end of the band shown in Figure [6] 
and conservatively apply it as a bound for all widths Tr/tiir < 0.25. 

3.4 Top pairs 

Searches for dijet resonances constrain the decays of R into the five lightest quarks. To 
constrain ■ BR(i? — > tt) independently, we use the dedicated CMS searches (56 57 for 
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Figure 8: Limit on g\ ■ BR(R -> 11) (solid) 
from the CMS search for dilepton reso- 
nances [49] and on g q ■ BR(i? — > tt) (dotted) 
from the corresponding ditau search [51]. 



Figure 9: Limit on g\ ■ BR(i? -> WW) 
from combined searches for Higgs decays to 



WW 52 



tt resonances. The resulting bounds are shown in Figure [7j 

We note that for family independent couplings and assuming itlr > 2m t , we can 
always use a bound on g q ■ BR(_R — > jj) to infer a bound on g q ■ BR(i? — » tt) using the 
relation 

BR(R tt) = V 2 l ' R ■ BR(R -> jj) . (10) 

A + 6 9dl9 u 

For (7 U ~ gd-, the bound on (7^ • BR(i? — > tt) inferred from the dijet limit is comparable 
to the direct bound from top pair searches. As expected, for g u <^ g d , the inferred 
bound on g 1 ■ BR(i? — > tt) is much stronger than the direct one, while for g u ^> g d , we 
can actually invert the equation above to obtain a bound on g 2 q ■ BR(i? — > jj) from the 
bounds on g 2 q ■ BR(_R — > tt). 

3.5 Dilepton resonances 



We consider the recent CMS search for dilepton resonances 49 with £ = (e, /1) as well 



as the search for ditau resonances 51 and show the resulting limits in Figure [8] If 
we assume family independent couplings, we can obtain a stronger bound on the ditau 
channel making use of the relation BR(rr) = BR,(££)/2. 
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100 150 200 300 500 700 1000 1500 
m R I GeV 

Figure 10: Limit on g\ ■ BR(R -»■ ZZ) 
from combined searches for Higgs decays to 
ZZ [52] (blue, solid) and from a dedicated 
search for high-mass ZZ resonances 
(green, dashed). 
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Figure 11: Limit on g\ ■ BR(R -> Z#) 
from the ATLAS search for associate pro- 
duction of Higgs and Z for a Higgs mass of 
m H = 125 GeV [55]. 



3.6 jyjy and ZZ 

Searches for WW and final state have been performed in the context of Higgs 



searches at the LHC 52 , 54 and as a dedicated search for high-mass resonances 53 



The results of these searches are shown in Figures [9] and [TOj We observe that the Higgs 
search limits from WW and ZZ on the R couplings are significant even for relatively 
large tur. The reason is that for a vector boson, the coupling involves a derivative, so 
it is enhanced compared to e.g. the SM Higgs coupling for large masses. 



3.7 ZH 

Recent results from the LHC exclude the SM Higgs with a mass between 130 and 

while the combined LEP2 results exclude 
In the allowed mass range, we 



550 GeV at the 99% confidence level 52,54 



63 



it below 114.5 GeV at the 95% confidence level 
have upper limits on the cross-section pp — > ZH from searches for Higgs production in 
association with a SM vector. We take the limit from pp — > ZH — > vvbb from 55 , 
giving <7 Z h < 1.7. . .2.0 pb depending on the value of the Higgs mass in the allowed 
range. We show the resulting limit for m# 



125 GeV in Figure 11 



12 



4 Implications for dark matter direct detection 



In this section, we apply the collider bounds obtained above to constrain DM-nucleon 
interactions mediated by R. We can calculate the direct detection cross-section by 
integrating out both R and Z to generate the corresponding effective operators. In the 
case of Dirac DM we obtain 

Cf = bjwxf-ff + &/X7m7 5 x/V7 5 / , (11) 

where we have neglected terms that vanish in the non-relativistic limit and defined the 
effective couplings 

A,V A,V A,V A,V 

bf Y = b% + b% = 9 -^p- + h*Ap_ . ( i2) 

Unless 6 V is very small compared to b A , the direct detection cross-section will be 
dominated by the effective vector-vector interaction between the DM particle and nu- 
cleons (p, n) given by 

4 = ftxi&pfp + fZxi.XnYn ; /* = 2b Y u + 6j , ft = 2b Y d + b Y u . (13) 
In the case of the complex scalar DM we have, similarly, 

£f = «/ J W . ^ = ft^PivP + f$Ji*w , (i4) 

where 

a] = a] R + a] z = g ^ + g -^- f} = 2$ + a? , f+ = 2$ + <% . (15) 

Because of the conservation of the vector current, there is no contribution of sea 
quarks or gluons to the effective couplings. For both Dirac and complex scalar DM we 
obtain the DM-nucleon cross-section 

ctn = I^nIn/ 71 ; where N — p, n . (16) 

In the following, we will consider two different possibilities for generating effective 
interactions of nucleons and DM particles. First we consider the case where R has 
sizeable direct couplings to quarks and all other couplings are arbitrary. Afterwards we 
consider the case where the interaction state X corresponding to the mass eigenstate 
R couples only to the DM particle and couplings to SM particles are generated only 
via mixing. An example would be the 'dark' Z', where R is the gauge boson of a new 
U(l) under which only the DM particle is charged. In this case, we can use collider 
bounds to directly constrain the mixing parameters, and therefore the direct detection 
cross-section. 
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4.1 Direct detection through direct couplings 

Let us start with the general case where R can have arbitrary couplings to SM particles. 
The only assumption we make is that R has a sizeable branching into quarks. This 
assumption is important for two reasons. First, we want to exclude the case where g q 
is so small that the total number of i?-particles produced at the LHC is insufficient to 
give a detectable monojet signal. We will come back to the case where g x 3> g^ R , gJ R 
m Section 0~1 

Second, this assumption ensures that DM direct detection is dominated by R- 
exchange, with Z-exchange giving only a negligible contribution. This case is interesting 
because it allows a ratio f n /fp significantly different from the one for Z exchange. In 
fact defining y = gi R jgi R we obtain f n / f p — (y + 2)/(2y + 1), which can in principle 
take any arbitrary value. 



We then obtain from Equation ( 16 ) 



„ (2y + lyA MM (g&) < (2 + 1)2 /4^ 

7r m R re m R 



For m x <C m R we can use Equation (24) to obtain 



a p < 12(2y + l) 2 x \ - gj • BR(i? inv 



m R 



As discussed in Section [3j monojet and monophoton searches at the LHC provide a limit 
on g\ ■ BR(_R — > inv), so that we obtain a bound on the direct detection cross-section if 
we can constrain T R . 

Of course, if we allow decays into new states that give complicated experimental 
signatures, we can make T R arbitrarily large. Therefore, we will now assume that all 
new states are either SM particles, or remain invisible, i.e. escape the detector without 
decaying into visible particles. In that case, we can combine the bounds from Section [3] 
to obtain an upper limit on and, assuming family independent couplings, constrain 



r R . The resulting bounds are shown in Figures 12 and 13 

With present data, we can constrain g^ and T R only in the range 300 GeV < m R < 
600 GeV. However, the only decay channel that is presently not available above 600 GeV 
is i? — ?- WW. We simply assume that upcoming searches for this decay mode will give 
bounds comparable to the current bounds for R — > ZZ. Consequently, we assume that 
the bound on g 2 q ■ BR(i? — > ZZ) also applies to g 2 q ■ BR(i? — > WW) so that we can 
extend our analysis up to 1200 GeV. Even a somewhat weaker bound on g 2 q ■ BR(i? — > 
WW) would not change our results dramatically, because decays into WW give only a 
subdominant contribution to the total width of R. 
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Figure 12: Limit on gj from the combination 
of all experimental bounds (see text). 




300 500 700 1000 1500 

m R I GeV 

Figure 13: Limit on Tr from the combination 
of all experimental bounds. 



We observe that up to itlr ~ 1000 GeV, Tr/itlr remains sufficiently small that the 
NWA stays valid, which is an important consistency requirement for our treatment. 
We can therefore use the limit on Yr from Figure 13 to calculate an upper bound on 
the direct detection cross-section. The resulting bounds on the direct detection cross- 
section both as a function of mediator mass and as a function of DM mass are shown 
in Figure 14 We observe that we can exclude a cross-section of a p = 2 x 10~ 40 cm 2 over 
the full mass range 300 GeV < rriR < 1000 GeV, as long as m x <C tur. 

So far, we have only considered standard spin-independent interactions of DM. Many 
other possibilities have been considered in the literature, e.g. spin-dependent interac- 
tions, momentum-dependent interactions, inelastic DM, and effective couplings with 
fn 7^ fp (see e.g. (64] ). Typically, these interactions strongly suppress scattering in the 
non-relativistic limit, while the results from LHC searches are not significantly affected. 
As a result, the LHC bounds will become much stronger compared to any exclusion 
limits (or claimed signals) from direct detection experiments. To illustrate this point, 
we show in Figure [15] the LHC bounds for spin-dependent interactions, as well as spin- 
independent interactions with f n /f p = —0.7. 

To conclude this section, we discuss how our results would change for complex scalar 
DM (e.g. a scalar technibaryon 36,37 65 66]). In fact, all the experimental limits can 
be applied in complete analogy, the only difference being that for identical couplings the 
partial width for decays into scalar DM is smaller by a factor of 4, cf. Appendix [A] As a 
consequence, the bounds on the direct detection cross-section will be weaker by a factor 
of 4, excluding a p > 8 x 10~ 40 cm 2 for m x < m R and 300 GeV < m R < 1000 GeV. 



15 



10 



-38 



10 



10 



-41 



-42 




300 



500 700 
m R I GeV 



DAMA 
CoGeNT 



fJf P = 
1000 



1(T 

io- 

10 
10 
10 
10 
10 
10 



39 



45 




CoGeNT 
DAMA 
XENON 100 
CDMS Gc low 
LHC combined bound for 
300 GeV < m R < 1000 GeV 
Contact operator bound 



m x I GeV 



Figure 14: Left: Bound from LHC data on the direct detection cross-section as a function 
of the mediator mass m# with the DM mass m x = 10 GeV. Right: Bound on the direct 
detection cross-section from LHC limits as a function of the DM mass m x compared to the 
results from various direct detection experiments. The width of the blue line corresponds to 
the change of the bound as the mediator mass is varied between 300 and 1000 GeV. For larger 
or smaller mediator masses, the bound would become weaker. As an example, we show the 
bound obtained from the contact operator if the mediator can be integrated out |18|. 



4.2 Direct detection through mixing 

Now we consider the case where R is the mass eigenstate that corresponds to the gauge 
boson X of a new U{l)x gauge group. X is then described by an effective Lagrangian, 
which includes kinetic mixing and mass mixing (see also Appendix [B]) 

C = Csm - \x^X, v + \m\X»X» - m x xx 



4 



1 



- - sine B^X^ + 6m 2 - f^X^X 



(19) 



We assume that the interaction eigenstate X couples only to the DM particle \ (as- 
sumed here to be a Dirac fermion) with strength /Y, and has no other couplings. In 
particular, we assume that other hidden sector states - even if present - give a negligible 
contribution to the total width of R. 

The mass eigenstate R then picks up SM couplings from mixing as described in 
Appendix [Bj At the same time, kinetic mixing introduces a coupling between \ and 
the Z-boson. After integrating out both R and Z, the resulting coupling constants for 
the effective interactions between DM and nucleons are, in terms of the Lagrangian 
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m R I GeV m R / GeV 

Figure 15: Left: Bound on the spin-dependent direct detection cross-section from LHC limits 
as a function of the mediator mass trr with the DM mass fixed to m x = 10 GeV. Right: Same 
but for spin-independent interactions with the isospin-violating couplings f n /fp = —0.7. 



parameters in Equation (19) 
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(20) 



where ^ is denned in Equation (39) and we abbreviated sin 8 = se, cos 9 = Ce,tan8 = tg. 
In the case where the mass mixing parameter Sm = 0, we obtain f n ~ 0, i.e. photon-like 
interactions. 

Because of mixing R can be produced directly in p-p collisions, so we can use LHC 
data to constrain the mixing parameters and therefore obtain bounds on the direct 
detection cross-section. In addition to the LHC bounds, we also have LEP bounds on 
the kinetic mixing parameters sine and £ and on g x z- In order to satisfy electroweak 



precision tests (EWPT) we must require that 23 

aS = 4^c^ v sw^e, and 

a T = ^ 2 (m^/m| - 1) + 2^s w t e , 



(21) 
(22) 



are within their experimental limits. Moreover, from measurements of the Z invisible 
width, we know that 



{g x z? < o.oos . 



(23) 



Note, however, that new physics might well give additional contributions to the 5* and 
T parameters, which can modify these bounds. 
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Figure 16: Bounds on the mixing parameter sine and the direct detection cross-section a p for 
fY = 0.1 and 5m = 0. The grey shaded region in the second plot corresponds to sine > 0.8. 
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Figure 17: Bounds on the direct detection cross-section for different choices of fY and 5m. 
The grey shaded regions correspond to sine > 0.8. 



For fixed fY and 5m we can calculate all coupling constants and therefore the partial 
decay widths and branching ratios of R as a function of tur and sin e. The bounds from 
Section [3] can then be interpreted as constraints on sine as a function of tur. These 
constraints directly correspond to limits on the DM direct detection cross-section. Our 



results are presented in Figure 16 



We observe that the bound from WW does not get stronger compared to the other 
bounds as we increase tur, despite the enhancement of the branching ratio of R into 
WW. The reason is that for small kinetic mixing, the coupling of R to quarks and 
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leptons is proportional to cos£t e ~ e, while the coupling of R to WW is proportional to 
sin£ rj esin^vF^|/ m jj- Thus, the partial width r picks up an additional factor of 
m\jm A R which precisely cancels the enhancement from the derivative interaction. 

In the case where /Jf = 0.1, the LHC gives strong constraints on the direct detection 
cross-section, comparable to the best bounds from current direct detection experiments. 
This conclusion does not change significantly, if we include a mass mixing term (see 
Figure 17). However, such a mass mixing will enhance £ compared to e, so the bound 
from WW becomes stronger compared to the bound from dileptons and monojets. 

Increasing /Y relaxes all bounds from the LHC since smaller quark couplings, and 
therefore smaller mixing parameters are required for the same direct detection cross- 
section. At the same time the invisible partial width of R is increased so that decays of 
R into SM particles are additionally suppressed. Still, even for = 1, we can exclude 
a direct detection cross-section above 10~ 41 cm 2 , because larger mixing parameters are 
excluded by both dilepton and monojet constraints as well as EWPT. 



5 Discussion 

In this paper, we have considered a new vector state R, as the dominant mediator 
of the scattering between DM and nuclei, relevant for direct detection experiments. 
We have demonstrated that the LHC can significantly constrain the scattering cross- 
sections over a wide range of mediator masses. In the case of fermionic DM and assuming 
that the vector mediator has direct couplings to quarks, we can exclude cross-sections 
o v > 2 x 10~ 40 cm 2 for 300 GeV < ttlr < 1000 GeV as long as m x <C ttir, with similar 
limits for scalar DM. When couplings to quarks are introduced only via mixing, the 
resulting bounds on the direct detection cross-section can be much stronger. However, 
it is possible to suppress the collider constraints by making the DM coupling very large 
compared to the quark couplings. Our approach allows a very general interpretation 
of the results from MET searches at the LHC, because we have not assumed that the 
mediator is heavy enough for an effective operator analysis. Nevertheless, we have made 
other assumptions in our analysis, which we discuss below. 

First, the DM particle has been assumed to have a mass m x <C w,r. When m x be- 
comes larger than R can no longer decay into DM particles and therefore monojet 
signals can only arise when the mediator is produced off-shell. Hence bounds from mono- 
jets become much weaker. We do not consider this case further, as LHC bounds on the 
interactions of heavy DM are significantly less stringent than those obtained from direct 
detection experiments. 

Second, we have assumed family independent fermion couplings (but we allow dif- 
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ferences between up-quark and down-quark couplings, and between charged lepton and 
neutrino couplings). Coupling the mediator predominantly to the third generation would 
slightly relax the bound from dileptons, but make the already quite strong bound on 
top-quarks even stronger. At the same time, doing so would strongly suppress the direct 
detection cross-section. Consequently, allowing family-dependent couplings would not 
significantly change our conclusions. 

Third, we have only considered the contribution from DY production of R. This is 
certainly justified in the case where couplings of R to standard model particles arise from 
kinetic mixing. For large ttir, however, the coupling of R to WW is not well constrained, 
so it could in general be possible to enhance the cross-section of VBF production. By 
neglecting additional contributions to the production of R we give of course a more 
conservative bound. 

Fourth, there is only a relatively limited range of mediator masses, viz. 300 - 
600 GeV, where we have bounds on all possible SM decay channels. Below 300 GeV, 
it will be difficult to constrain some of the couplings of R using LHC data, mostly be- 
cause of QCD backgrounds. Instead, EWPT will become more important (see [33]). To 
extend the range of mediator masses above 600 GeV, we have simply assumed that the 
bound on BR(_R ->■ WW) above 600 GeV is identical to that on BR(R ->■ ZZ). Since 
the total width of the mediator is not very sensitive to this bound, we expect our result 
to be only mildly affected if the observed bound is significantly weaker. 

Finally, and most importantly, we have assumed that we can treat R as a narrow 
resonance and we make extensive use of the NWA to separate the production of R 
from the subsequent decays. This factorisation is not valid for broad resonances, where 



we expect a significant contribution from off-shell mediators (see 45 for a discussion). 
Moreover, most bounds only apply for a narrow resonance. For example, limits for dijet 
resonances have only been published for peaks with Y/M < 0.3. 

We observe that for tur < 1000 GeV, present experimental results give the constraint 
^ n/ m R < 0.25. For these values, we expect the NWA to be accurate within a few 
percent. As tur increases, the bounds on the individual coupling constants, and therefore 
the bound on Yr, become weaker. For tur > 1 TeV, couplings can be of order unity for 
which the width becomes so large that neither experimental limits nor the formalism 
presented in this paper can be applied. With increasing luminosity at the LHC, we 
expect to be able to extend our treatment to larger mediator masses, while at the same 
time obtaining stronger bounds on the direct detection cross-section. 

For dark matter masses m x <C uir, dark matter annihilation will only involve off- 
shell mediators. If the mediator couples mostly to quarks and dark matter, it is therefore 
possible to describe annihilation in terms of the same effective operators as direct de- 
tection. Consequently, one can directly translate between direct detection cross-sections 
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and annihilation cross sections (see for example 67,68]). In the case of a vector mediator 
and light dark matter, direct detection cross-sections of 10~ 40 cm 2 or smaller correspond 



to an annihilation cross section significantly below the thermal cross section 14 , 69 



The LHC bounds on the decay channels of R therefore imply that the annihilation 
cross section of dark matter into quarks is well below the present experimental sensitivity 
for dark matter indirect detection in all relevant channels including the antiproton flux 
from PAMELA [TO] , diffuse gamma rays from FERMI-LAT 71 73 and neutrinos from 
. 3 Even more restrictive experimental bounds in the future can be evaded 



the sun 74 



either if dark matter is asymmetric or by appealing to additional uncertainties as in 76 
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To obtain the required dark matter relic density in our framework, we must assume 
either the presence of additional mediators in the early universe in order to avoid over- 
production of DM or additional couplings of R to new hidden sector states. In fact, 
coupling R to new hidden sector states that decay into SM particles with more compli- 
cated experimental signatures is an interesting possibility to evade experimental limits 
on T^j. Thus there are good reasons to carry out experimental searches 
example) for such states. 



see 



78 for an 
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'However, if R has only small coupling to quarks, there could be a monoenergetic gamma ray signal 



from annihilation of dark matter into Z gamma as in 75 
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where xz = (rnz/m R ) 2 , xr- = (m# /rn R ) 2 , and \(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. 
Note that in the latter formula g§ H has mass dimension one, following our conventions 
in the main text, and the expression agrees with that given in e.g. [79] . 
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B Coupling structure from mixing 



In this appendix, we discuss how the mass eigenstate R arises from the mixing of an 
interaction eigenstate vector X with the SM U(1)y B field and the neutral component 
W 3 of S77(2)l weak fields. We first consider the most general case and calculate the 
effective coupling constants defined in Section [2] in terms of the fundamental couplings 
and the entries of the mixing matrix. We then calculate the mixing matrix that arises 
from kinetic mixing and mass mixing of gauge bosons. 

Following the notation in [33] we write the general mixing matrix as 

4\ / N u N 12 N 13 

1 = 1 N 21 N 22 N 23 | ( Z^ | . (31) 



X„ / V N 3l N 32 Ni 



Here A, Z are the physical photon and neutral massive gauge boson fields of the SM. 
The couplings of R to SM fermions are given in terms of the mass mixing matrix, as 4 

9uR = ~WN a + 3gN 23 ) - f^N 33 , 9uR = \(g'N 13 - gN 23 ) - f^N 33 , 

9m = ^(g'Nis + SgN 23 ) - f^N 33 , g m = -\{g' N l3 - gN 23 ) - /j^s , 

9eR = \(3gNi 3 + gN 23 ) - / e v iV33 , <4 = -\{g'N 13 - gN 23 ) - f*N 33 , 

9uR = \(9'N 13 - gN 23 ) - fiN 33 , tfi R = -\{g ! N l3 - gN 23 ) - f^N 33 , (32) 

where the numerical coefficients are determined from the hypercharge and weak quantum 
numbers of the SM fermions and / V,A denote the direct couplings of X. Similarly, the 
effective vector and axial couplings of R to the DM particle are given by 

9xR = fx N ^ > tin = f£N 33 or 9<j>R =UN 33 (33) 

depending on whether the DM particle is a fermion or a scalar. 



Finally, the couplings of R to SM bosons are given by 28 



9wwi = 9ww2 = 9N 23 , Qzwwx = ~9 2 N 22 N 23 , g\ wwx = -g 2 N 21 N 23 , (34) 
9zh = \{g'Nn - gN 22 ){g'N l3 - gN 23 ) , g§ HH = -(g'N 12 - gN 22 ){g'N 13 - g)N 23 ). 

(35) 



4 We use g and g' to denote the fundamental gauge couplings of SU(2)l and U(l)y, which will be 
different from the observed ones, g and g' . 
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Since P-violating couplings of gauge bosons are absent in the SM, the corresponding 
couplings of R cannot be introduced by mixing alone. 

We now assume that X is the gauge boson of a new U{l)x gauge group and follow 
the discussion in [33] of an effective Lagrangian which includes kinetic mixing and mass 
mixing (see also |23|) 



C — Csm — -X^X 



4 



-m%.X u X' t 
2 



m x xx 



- \ sine + 5m 2 Z,X» - ^ fJX^f^f - f^xi.X 



(36) 



As in [38], we assume that the U(l)x is broken by an additional Higgs field and X 
acquires the mass m^. We will not discuss the implications of this additional Higgs 
field and its potential mixing further. We define Z = c^W 3 — sw-B, where % (cw) is 
the sine (cosine) of the (fundamental) Weinberg angle. 

The diagonalisation of the above Lagrangian is discussed in detail in e.g. [23]. The 
field strengths are diagonalised and canonically normalised via the following two con- 
secutive transformations 





where 







— 2c £ (Sm 2 + m|sws e ) 



.4, 



R,. 



%c 2 + m%s 2 v s 2 + 2 5m 2 s w s e 



m 



x 



in 



(37) 



(39) 



Multiplying the two matrices, we obtain the coefficients iV™ so that we can calculate 



the couplings of R using Equations (32 35). As discussed in 28 33 , the fundamental 



parameters and % are constrained by the requirement that the physical Z mass 
and the Weinberg angle come out in accord with experiment. 

At colliders we are directly sensitive to the couplings g U} d, which determine the 



production cross-section of R (see Section 2.2 ). We show these in Figure 18 as a function 



of mn for different values of the kinetic mixing parameter e and different values of the 
mass mixing parameter 8m. 
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Figure 18: The couplings g\ (solid lines) and g 2 d (dashed lines) as a function of trr for e = 0.01 
(green), e = 0.1 (blue) and e = 1 (purple). In the left plot, interactions are induced by kinetic 



mixing only via the Lagrangian given in Equation ( 36 ) , while in the right plot we have included 
a mass mixing of 5m = mz/2. 




sin e sin e 

Figure 19: The effective couplings f p (solid lines) and f n (dashed lines) obtained from the 
Lagrangian given in Equation (36) as a function of sine for mR = 300 GeV (green), = 
600 GeV (blue) and mR = 1200 GeV (purple). In the left plot, interactions are induced by 
kinetic mixing only, so f n = 0, while in the right plot we have included a mass mixing of 
5m = mz/2- 



If there are no direct couplings to quarks, we can calculate the effective DM-nucleon 



couplings f p and /„ in terms of e, 5m and f^, cf. Equation (|20|). We show these couplings 
as a function of sine for different values of 5m in Figure 
obtain f n /f p = 0, while for 5m ^ 0, f n /f p < 0. 



9 Note that for 5m = 0, we 
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